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Antipodally  symmetric  distributions  for  orientation  statistics 


Michael  J.  Prentice 
Department  of  Statistics 
Edinburgh  University 
Scotland,  U.K. 


SUMMARY 

The  conventional  antipodally  symmetric  Bingham  matrix  distribution 
on  the  Stiefel  manifold  is  generalised.  Large  sample  maximum  likelihood 
estimation  and  uniformity  tests  are  discussed,  and  a  parametric  model 
for  axial  orientations  (X-shapes)  is  suggested.  A  generalisation  of  the 
Khatri-Mardia  matrix  distribution  is  developed  to  provide  a  model  suitable 
for  hybrids  (T-shapes).  Beran’s  results  on  exponential  models  for  direc¬ 
tional  data  are  extended  to  orientation  statistics  to  provide  regression 
estimators  and  goodness-of-fit  tests  as  alternatives  to  maximum  likeli¬ 
hood  estimation  and  likelihood  ratio  tests. 
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1.  Introduction 

Following  Downs  (1972) ,  wa  define  an  orientation  statistic  as  a 
rigid  m-frame  in  R^  (m^p)  i.e.  an  m*p  matrix  X  s.t.  XX'  *C  , 
idiere  C  is  an  m*m  symmetric  positive  definite  matrix  specifying  the 
angles  between  the  rows  of  X  .  Without  loss  of  generality  we  suppose 

C»I  since  all  methodology  far  C  - 1  can  be  extended  trivially  for 

m  ■ 

general  C  .  We  define  an  L-shape  to  be  a  rigid  m-frame  of  signed 
directions  (a  conventional  orientation  statistic) ,  an  X- shape  to  be  a 
rigid  n>- frame  of  axes ,  and  a  T- shape  to  be  a  hybrid,  an  m-frame  of  m^ 
axes  and  ■*  m  -  signed  directions  in  R*  .  All  methodology  so  far 

published  is  suitable  for  L-shapes.  We  assume  that  p>,3  since  when 
p  ■  2  every  type  of  orientation  statistic  reduces  to  either  a  direction 
or  an  axis  an  the  circle. 

The  van  Mises-Fisher  matrix  distribution  provides  a  suitable  unimodal 
pdf.  for  orientation  statistics.  Maximus  likelihood  estimates  and  likelihood 
ratio  tests  have  been  developed  by  Downs  (1972) ,  Xhatri  and  Hardia  (1977) 
and  Jupp  and  Hardia  (1979).  The  conventional  fully  parameterlsed  Bingham 
matrix  distribution  (Xhatri  and  Hardia,  1977,  (7.2)  with  6  -O  )  is  the 
obvious  analogue  on  the  Itiefel  manifold  of  Bingham's  antipodally  syanetric 
distribution  on  the  sphere  (Bingham,  1974) .  Maximum  likelihood  estimation 
and  tests  of  randomness  against  a  special  case  of  this  distribution  have 
been  treated  by  Jupp  and  Mardla  (1979)  and  Hardia  and  Xhatri  (1977) .  In 
Section  2  we  state  the  corresponding  MU  results  Jot  the  fully  parameterlsed 
2*  -modal  Bingham  matrix  distribution.  We  extend  these  results, -to  a 
generalisation  of  the  Bingham  matrix  distribution  in  Section  3,  and  obtai|n 

i 

a  new  'Bingham'  statistic  as  a  large  sample  test  of  uniformity  against. 

» 

s 

general  antipodally  symaetric  alternatives,  a  large  sample  approximation 


to  the  likelihood  ratio  statistic.  In  Section  4  s  special  csss  of  this 
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distribution  is  suggested  ss  a  suitable  nodal  for  X- shapes,  and  an  axial 
Binghaa  statistic  is  obtainsd  for  a  tast  of  uniformity.  Ths  generalised 
Xhatrl-t!ardia  Matrix  distribution  is  usad  in  S action  5  to  provide  a 
suitable  par ana trie  nodal  for  T-shapes,  and  a  hybrid  Rayleigh-lingham 
statistic  is  obtained  as  a  large  sample  tast  of  uniformity.  Beran’s 
results  on  rotationally  invariant  exponential  aodals  for  directional 
data  are  axtandad  to  the  8  tie  fa  1  nanlfold  in  Section  6*  to  provide 
regression  estlaators  and  goodness  of  fit  tests  within  t be  generalised 
Xhatri-Hardia  family  of  distributions#  as  alternatives  to  the 
aoaputatlanelly  inconvenient  naxinun  likelihood  estlaators  and  likelihood 
ratio  tests. 

2.  the  BioUhaa  astrix  distribution 

Zf  X  is  an  amp  astrix  randan  variable  (a£p)  with  pdf 

C2»)“^  |x|*  M*  etr(-«*(X-u)V(X-w)')dX  (2.1) 

dart  X  is  wm  aynetrlc  positive  definite#  and  V  is  pg  syanetrlc 
positive  definite#  then  if  ii«Q^  #  an  sap  astrix  of  seres#  the  pdf 
of  X  #  conditional  on  XX'  -Im  is  (Xhatrl  and  Hardla#  1977#  (7.2)) 

etrC-ym-1  I  -J#XVX')  CdX3  (2.2) 

m 

where  (dX)  denotes  the  unifora  distribution  on  the  8 tie f el  nanlfold 
0(a#p)  #  and  Y  ■  Y  (X  •  V)  is  a  no  realising  constant  which  depends  only 
upon  the  diagonal  aatrlees  0^#  of  eigenvalues  of  X  and  V  . 

t  Y 

Series  expansions  for  a  #  a  hype  rge  one  trie  function  of  two  astrix 
servants#  have  bean  given  by  Srivastava  and  Carter  (1980).  For  reasons 
which  will  bacons  apparent#  we  denote  the  distribution  (2.2)  >(X#a#p,KSV) . 
Xts  pareaeter  space  has  diaension  (*^1)  a  (pJl)  -2  when  l«,n£p-l  #  since 
B(X#a#p#  («x)  •  (a  1  gaSlp))  ■  B(X#a#p#xgg)  for  all  real  scalars  ado#B 
(cf.  Binghaa#  1974#  Zanaa  2.1).  Biailsrly#  when  s»p 
8(X»P#P»  •  (•  1  ■  B (x#p#p#  xiv)  ,  where  without  loss  of 


4 


generality  ve  nay  assume  trace  (X-V)  ■  O  ,  so  that  the  parar«ter  apace  has 

dimension  2(^1)  -  3  .  Where  convenient  to  identify  pa  masters  uniquely 

vs  shall  assuas  that  w  *0  *pps°*  »"P)f  °*r  if  dealing 

in  spectral  decompositions ,  that  the  entries  in  «©K  and  Sty  are 

in  decreasing  order  with  |Il),.»0  (and  (D_)  *0  if  m*p  ).  He 

v  11 

issues  that  tbs  elscents  of  Dj,  and  Dy  are  distinct  so  that  K.  has  a 
unique  natrix  Q  cO(i*  ,a)  of  eigenvectors  and  V  has  a  corresponding 
matrix  MeO(p-l,  p)  .  Proa  Theobald  (1975,  Theorea  1),  the  distribution 
(2.2)  has  2s  nodes  at  the  points  X»Q'M^  ,  where  is  the  a'  xp 
astrix  of  the  first  a*  >aln(a,p-l)  ooluans  of  M.  The  aultlplicity 
2  arises  from  the  possible  sign  changes  of  the  columns  of  Q  and  . 

Z  aa  grateful  to  Dr  Theobald  for  first  drawing  ray  attention  to  this. 

A  random  saaple  Iji  . * . ,  on  0(a,p)  from  the  distribution 

(2.2)  has  log  likelihood 

-ay  -  *m  trace (K#v)t  (2.3) 

-1  " 

there  l*a  £  X^^fX'  •  (y^|q)  ,  1<X*1<P  .  The 

distinct  eleaanta  of  T  ,  excluding  thoee  for  whldi  (k,l)  »  (p,p) 

(and  also  thoee  far  which  (j,q)  •  (p,p)  when  n«p)  are  sufficient  but 
not  minimally  sufficient  unless  a  ■  1  (the  case  of  directional  data) , 
>«.  .  <■•*>  -  (*;l> .  rj1)  *  -  a  ,  with  equality  only 

whaa  «•!  ,  and  v(p,p)  -  t***1)  -  2  (*JX)  ♦  l^^1)  -3  .  The  following 
is  a  consequence  of  Berk  (1972) ,  and  elementary  calculus. 

Theorea  2.1 

(a)  Aar  sufficiently  large  n  ,  there  exist  unique  MUXs  K,V  subject  to 

the  condition  V»©  (and  fc  *0  and  trees  (X~V)  *0  whan  ■•»)  of  the 
PP  PP 

parameters  in  l(X,i,p,XI7).  The  MU  are  the  eolations  of 
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.  a  _1  Q 

where  W.(V)  -  n  A  I  X4VX*  and  *,(X)  -  n  I  X!  XX  . 
1  i-i  4  1  2  i-1 


(b)  (Spectral  nrslcn)  Zf  K-Q'D^fl  and  V-M'iyi  are  unique 
spectral  decompositions  ef  K  and  V  #  where  Q  cO(n*  ,a) ,  Me  0(p-l,p) 
end  D^,  Dy  are  respectively  a*  *  a'  and  (p-1)  *  (p-1)  diagonal 
aa trices  (and  trace  (D^-Dy)  -0  if  n-p)  ,  then  tor  sufficiently 
large  n  with  probability  1  there  exist  unique  MLEs  ft#  S^#  fiy  of 


Q,  n,  Oj^«  Qy  given  by  the  unique  spectral  dscoaposltlons 

(V)  -Q'Dy  Q,  Wjdb  -M'fiy  ft  and  the  equations 


(*) 

(*) 


<  W  -  <v  V 

«v  v  -  lvV 


l 

-%  ^  . 
^2 


Me  offer  no  algorithm  for  the  evaluation  of  ft  and  ft  •  Given  suitable 
Initial  approximations  X^,  ftQ  it  should  be  possible  to  construct  an 
Iterative  procedure,  given  tractable  series  e^ansiens  for  y  aad  Its 
first  derivatives  (see  Binghaa,  1976)*  but  considerable  oonputatloaal 
effort  Is  necessary. 

App  and  Nardia  (1979#  Theorems  1(b)  aad  4)  have  obtained  the 
analogue  of  Theorem. 2.1  for  the  special  case  KXfSrPtZ^PV)  which  has 
parameter  space  of  dimension  %(p-l)(p*2)  provided  m  <p  .  we  note 
that  B(X#p#p«ZpgV)  Is  the  uniform  distribution  so  their  results  are 
invalid  for  equere  orlentatlona.  The  log  likelihood  of  a  saaple  fron 


(2.4) 


is  -n  y(Ia§V)  -ta  trace  VT* 

»i  n 

where  T*  -  n  z  X.'X.  is  pap  aynaetric,  positive  definite  with 
i-1  1  1 

probability  1.  The  distinct  elements  of  Y*  ,  excluding  y^  say,  are 
ainlaally  sufficient  for  V  ,  subject  to  vppm°  **y.  A  einple  large 
sample  test  of  unifocal ey  against  alternatives  (2.4) ,  asyaptotically 
equivalent  to  the  likelihood  ratio  statistic,  nay  be  obtained  by 
generalising  Bingbaa's  Theoraa  5.2  (1974,  p.1206) .  He  obtain 


Yheoren  2.2 

On  the  null  hypothesis  of  uniforaity  the  statistic 
•B«  (traced*2)  -n*/p)  np(p~l)  (p +2)/2a(p-a)  is  asyaptotically  distributed 
as  x2  an  v(l,p)  -h(p-l)  (p+2)  degrees  of  flreedon. 


Proof  (A  simpler  version  of  Theaoren  3.2  below). 

Xf  is  the  -vector  of  the  distinct  elenents  of  T*  ,  then 
X*,  the  esyaptotic  varlanoe  aatrix  (Mar die  and  Xhatri,  1977  p.469)  of 
»V  has  rank  v(l,p)  and  a  generalised  inverse  (p(p-l)  (p+2)/2a(p-n)  ] 


-1 


bloekdiag  (lp»  2X^  )  .  Since  Y*  has  null  expectation  op  *  xp  ,  the 
result  follows  laws diately  (see  also  xhatri  and  Mardia,  1977,  p.471,  for 


alternative  derivation). 


*99£*£  1- 
2. 


S  is  undefined  vben  a»p 
u 

B(x,a,p,X«Xp),  the  other  obvious  special  ease  of  (2.2),  is 
the  uni  fora  distribution  for  all  ■  ,  l<«<p  . 


A  qmTAl  1  iiti  mi 

Consider  the  fully  paranetsrlsed  aultiveriate  noraal  nodal 
(2*)“^  |a|^  etr(-%  A(x-e)  •  (x-*)')dx 
and  the  oonditienal  distribution 


(3 


(3.2) 
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etrf-ytap)"1  1^,  -  h  AXgX’)(dX) 

obtained  when  XX'  -  I  and  u«0,_  .  Bara  A  is  apxop  synmetric 
poaltiva  definite  with  distinct  eigenvalues  a^x ...  »a  ,  and  y  my  (A) 
la  a  normalising  constant.  The  distribution  (3.2)*  denoted  B(X,m,p»A) 
is  antipodally  symmetric  and  has  paremetar  space  of  dimension  v(m,p)  , 
since  B(X,m,p,A)  -B(X,a,p,A  +9^  •  Ip)  for  all  real  mum  symmetric  matrices 
A^  ,  and  B(X,p,p,A)  -B(X,p,p,A+A1BXp*XpBA2)  for  all  real  p*p  syxssetric 
matrices  A^,  Aj  *  where  without  loss  of  generality  we  may  assume  trace  (A^  -  A, 
Where  convenient  to  identify  parameters  uniquely  we  shall  assume  that 
A  •  (•  jk,  £q)  «£■#  l£h,l£p»  satisfies  the  conditions 

a.,  _ "O  for  all  J,q,  and,  if  m«p,  a^  .“O  for  all  k,i  .  C 
3P,pq  **»*■? 

A  random  smople  from  the  distribution  (3.2)  has  log  likelihood 

-nr (A)  -ha.  trace  AX  ( 

with  T  -  (yikflq>  u  in  (2.3).  The  distinct  elements  of  T  ,  excluding 
thoee  corresponding  to  (3.3) ,  are  minimally  sufficient  for  A  .  From 
Berk  (1972) ,  and  differentiation  we  obtain. 


Theorem  3.1 


For  sufficiently  large  a  ,  with  probability  1  there  exists  a  unique 
MX  A  of  A  subject  to  (3.3),  which  is  the  solution  of 

-  -%*  . 

(M) 


As  with  theorem  2.1,  there  are  considerable  computational  difficulties 
associated  with  the  search  for  A  .  Xn  particular,  no  explicit  fora  of 
Y (A)  is  currently  available. 

A  Bingham  statistic  for  testing  uniformity  against  all  antipodally 
sy—  trie  alternatives  of  the  form  B(X,a,p,A)  may  be  obtained  as  follows. 


,)  *o  . 


.3) 


.4) 
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Theorem  3.2 

*•»  •“*  *'WW  •  0,1  th* 

null  hypothesis  of  uniformity  the  statistic 

Bnjp»n(p-l)Cl5(p+2)  trace  (Z^Z^)  -  trace  (Z^Zg)!  is  asymptotically 
distributed  as  x2  on  v(m,p)  degrees  of  freedom. 


Proof 


-1 


On  the  null  hypothesis  ?  has  expectation  p  1^  .  Consider  first 
the  case  a<p.  the  distinct  elements  of  Y  nay  be  written  as  a 

sector 

X- (<  XU> . <  W ' <^12> . <Vu>)' 

where  each  <Yjj>»  l£3  <.»*  i*  a  C1^1) -vector 

(yjl#lj#  y32*2j'  yjj,jj‘p  l'  ***'  y3p»pj '  yjl»2j'  yjl,33r  yJp-l#p3) 
and  each  <YjJl  *  •  i£J  * *l*i  is  a  p2 -vector 

tyjl*li*  yjp,pl'  yjl,2i'  y32,li'  *•*'  y3p-l,plr  y3PrP“l*>'*  °*ing 
Anderson  and  Stephens  (1972*  p.616)  It  follows  that  n*y  has  covariance 

matrix  I  ■  [ (p-1)  (p+2J  1  blockdiact (20^ Bblockdlao  (0^*21  ))» 

-1 

X_  •  blockdiaq (C  ■  »p  X  n  0  0)) 

Q  W  $ 

where  C  ■  X^  -  p  1  1— *  1^  is  an  e*m  matrix  of  ones*  and  D 


•fc  ;']■ 


Sines  C-1  “C 


follows  that 


and  C  ■  X  is  a  generalised  inverse  of  CM  *  it 
PP  P  PP 


l"- blockdiag ( (*i {p-1 )  (p+2)C  ^BblockdiagtX.lI  )) , 

■*»-p  P  (Pj 

X.  •  blockdiag{(p-l)  (p+2)  X  -X  M  BD.)) 

<j  p  $  1 

[(P*-l)r  (P- 1)1 
(P"l)  <p*-l)  *  ** 


Z  has  rank 


and  the  quadratic  form  a  y'l  y  reduces  to  B^  as  above*  it 
follows  that  B^  is  asysptotioally  distributsd  as  Xs  on  v(m*p) 
of 
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11m  only  ehano*  necessary  when  ■  ■  p  is  that  C  should  ba 

®#*“P 

replaced  by  X_  (■  C _ ) .  We  obtain  B _ as  above,  where  now  £  has  rank 

P  PP  PP 

v(p,p). 


Remarks  1.  Whan  p>_4  ,  B  is  pradsaly  the  'Bingham*  statistic 
obtained  by  Prentice  (1981,  (3.4))  from  consideration  of  the  (0,2) -th 
and  (0,l,l)-th  characters  of  th*  lrraduelbla  continuous  raprasaatations  of 
th*  rotation  group  0*(p). 

2.  whan  ■  ■  1  »  is  the  conventional  Bingham  statistic  for 

dir  actional  data,  since  than  Z„  ■  0  and  Z.  is  symmetric. 

z  -np  J. 


4*  Axial  orientation  statistics 

The  distribution  (3.2)  (or  (2.2))  may  ba  specialised  to  give  a  parametric 
modal  suitable  for  X- shapes.  We  require  a  probability  density  on  0(a,p) 
invariant  under  sign  changes  of  any  row  of  the  random  variable  X  .  this 
is  achieved  if  a  satisfies  the  conditions  •jjnq"0  **  j  +<l  •  We 
write  (X,m,p,E)  for  the  density 


V'iW 


(4.1) 


where  Y  "Y  (*)  Is  a  normalising  constant,  I  is  an  aspxp  array  with 

3th  layer  ,  1  <.  j<.n  ,  *  !*»*£?»  «4  Sj  1*  the  J-th  sow  of  the 

random  variable  x.  Since  B(mc)  (x,m,p,B)  -B(“t)  (X,m,p,I-»C*)  for  all  real 


aepxp  arrays  with  *1 


•3JU-3  • 


that  the  distribution  (4.1)  has 


:  an  layer  only,  it  foil 
of  dimension 


n(m*p)  «m  (*£*)  - n->i  m(p-l)  (p+2) ,  when  *«p  .  Xf  n»p»  ej,^  nay  lx 
the  more  general  form  e^+  f^  where  F»  (f^)  is  p*p  symmetric,  and 


be  of 


without  loss  of  gansrsllty,  f  -o  . 

PP 


for  square  axial  oriaatations 


of  the  distribution  (4.1)  has  iM mansion 


fl(p,p)  -p^1)  -  P^1)  -p  +  l-^(p-l)a(p4-2)  •n(p-l,p)  .  where  oonvenisnt 


to  identify 


uniquely,  we  shell  asa 
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•j  *0  for  all  j  ,  l<.  j£m,  and  If  m»p,  8^-0  alao.  (4.2) 

A  random  sasple  from  the  distribution  (4.1)  has  log  likelihood 

a 

-irr(E)  -*jn  trace  Z  E.V.  (4.3) 

j-1  3  3 

-1  “  (il  (11  •  (11 

where  Y  *  n  £  1  x.;  ,  x!'  representing  the  jth  row  of  X.  . 

The  distinct  elements  of  Y, ,  . . . ,  Y„  ,  excluding  those  corresponding  to 

l  n 

(4.2) ,  are  minimally  sufficient  for  E  .  As  in  Section  3  we  obtain 


Theorem  4.1 

For  sufficiently  large  n  ,  with  probability  1  there  exists  a  unique 
MLE  E  of  E  ,  subject  to  (4.2) ,  which  is  the  solution  of 


W 


a  *  "  %(E.»  »••»  Y  ,)  • 


(EsE) 


An  axial  Bingham  statistic,  suitable  far  testing  uniformity  against 
alternatives  (4.1)  aay  be  obtained  from  a  simplified  version  of  Theorem 
3.2.  From  consideration  of  the  asymptotic  null  distribution  of  the 
first  elements  of  n^y  ws  obtain 

Theorem  4.2 

On  the  null  hypothesis  of  unifomity  the  statistic 

-  ‘jn(p-l)  (p+2)  (trace (  Z  Y  2)  -a/p)  is  asymptotically  distributed 
V  j-1  3 

as  x2  on  n(m,p)  degrees  of  freedom. 


Remark  8  (Theorem  2.2)  -Hnp(p-l)  (p42)  (trace ((  Z  YJ*) -m2/p)/e(p-m)  , 

J-1  3 

provided  e <p  . 


-11- 


Bybrida 

A  parametric  nodal  suitable  for  T-shapes  nay  be  obtained  £ron  the 


generalised  Xhatri-Mardla  distribution 


-1 


Imp  *  SAX*X'  +  Ay#X*)  , 


etr(-1f  (up) 

from  (3,1),  conditional  on  XX*  ■  I  ,  when  m*0, 

2D  *  *1 


(5.1) 

We  require  that  the 


pdf  should  be  invariant  under  sign  changes  of  any  of  the  first  m^, 
rows  (the  axes)  of  X  .  This  can  be  achieved  by  requiring  that 


<n 


y  •  '•here  !•  an  n^xp  matrix  of  naans,  O  <  n^  ■■ 

and  that  A  satisfies  the  condition  a4b  .  ■  0  if  }  y<q  and 

min(j,q)  •  We  write  (X,a^,n2,p,u2,E,A2)  for  the  density 

“l 

exp (-Y “ **  trace (  £  ♦  trace  (Ajjij  SIX*  2J  A2X(2)  •X(2),) 


(5.2) 


where  y«y(E,A2)  is  a  nornalising  constant,  E  is  as  in  Section  4,  but 
with  only  n^  layers,  is  the  a^xa^p  submatrix  of  A  corresponding 

to  X^j ,  the  last  Wj  rows  of  the  random  variable  X  .  Since 

jOty)  (x,n2,n2,p,u2»E,A2)  (X,n^,m2,p,U2,E  +  E*,  Aj+AJ#*^) 

far  all  real  B*  as  in  Section  4  (but  with  only  m^  layers) ,  and  all 

real  symmetric  matrices  AJ  ,  it  follows  that  when  n-n1+n2<p, 

the  distribution  (5.2)  has  parameter  space  of  dimension 
.  mjP+l  n-+l 

♦  (n^m^p)  •n1(p^J')  ♦  (  2  )  -a^-  (  2  )+  Sjp  •  ritm^p)  ♦  vOa^p)  +WLp  , 


When  m«p  ,  the  more  general  result 

l0**  (X^,^,P,V2,E,A2)  -  {Xfmltm2tp,V2fZ-*t*tA2*  (AJ»Ip)  ♦  (X^»»)) 

obtains,  where  F  and  B*  are  of  the  more  general  form  in  Section  4,  and 
so,  the  distribution  (5.2)  has  parameter  space  of  dimension 

fta^m^p)  -  n(n£,p)  ♦v(n2,p)  ♦njp  ,  whore  m£  ■  nintn^p-Oj-l) . 

Where  convenient  to  identify  parameters,  we  shall  assume  that  Aj 
satisfies  (3.3),  and  that 

•jgp  ■  O  for  all  3,  11,3  im^,  sod  if  m^<(^«p,  1*0  (5.3) 
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A  randan  sample  from  the  distribution  (5.2)  has  log  likelihood 

"l 

-rry  (e,a2)  -  4n  trace  (  Z  +  n  trace (A2v2  •  X'^  -  AjT^j)  (5.4) 

where  Xj2j  in  the  n^xp  matrix  of  means  of  the  last  m2  rows  of  the 
data  matrices  X^,  . ..,  XQ  ,  and  Y  q)  is  the  m2^>  x m2^3  *ubinatrix  of  Y 
corresponding  to  the  last  m2  rows.  The  distinct  elements  of 
5c v  Y^,  ...  Y (2)  *  •xcludin9  those  corresponding  to  (3.3)  and 

(5.3),  are  minimally  sufficient  for  y2,  E  and  A2  . 

Theorem  5.1 

For  sufficiently  large  n  ,  with  probability  1  there  exists  a 
MBS  0  ■  (u2,A2,E)  of  0  ■  (u2,A2,E) ,  subject  to  (3.3)  and  (5.3),  which 


is  the  solution  of 


(*L  . . v 


A  large  saqple  test  of  uniformity  of  T-shapes  against  alternatives 
(5.2)  may  be  obtained  by  slight  modification  of  Theorems  (3.2)  and  (4.2), 
We  obtain  a  hybrid  Bay lei gh-Blngham  statistic 


On  the  null  hypothesis  of  uniformity,  the  statistic 


*t““*°*  yi)'*1"1  “  1”  "-O'- 

Mm  n  as  in  Theorsm  3.2,  is  asy^totlcally  distributed  as  x2  on 
+(m£,m2,p)  degrees  of  freedom. 


(ax) 

r  _  as 


in  Theorsm  4.2,  and 
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Proof 


*2  has  expectation 

V  ***'  ^  “d  *(2)  * 

quadratic  fora  in  BjP  +  ( 
matrix  of  rank  +  (a^,m2,p) 


Qb  and  is  asymptotically  uncorrelated  with 

Using  the  method  of  Theorem  3.2  we  obtain  a 
®2P+1  P+1. 

2  )  ♦  r2  )  variables  with  covariance 
,  which  reduces  to  as  stated. 

D.D-P 


Exponential  models  and  regression  estimators 

It  is  instructive  to  consider  these  matrix  Bingham  and  Khatri-ilardia 
distributions  within  the  context  of  the  rotationally  invariant  exponential 
family  of  distributions  for  orientation  statistics,  obtained  by  generalising 
Beran'a  (1979)  exponential  family  for  directional  data.  Provided  good 
multivariate  density  estimates  are  available,  the  obvious  analogues  of 
Beran'a  regression  estimators  ((1.10)  ibid.)  and  goodness  of  fit  tests 
(Section  5,  ibid.)  should  be  considerably  more  convenient  computational ly 
than  exact  HLEs  and  likelihood  ratio  tests. 

Consider  first  the  case  m  ■  p  ,  and  the  general  exponential  model 
exp(h(X)  -  y  (h) ) ,  h  c  •  ,  where  is  associated  with  the 

g-th  character  of  the  irreducible  continuous  representations 

of  0+(p)  ,  as  in  Prentice  (1981).  The  von  Mlses-Fisher 

matrix  distribution  spans  ^  ,  of  dimension  p2  ,  end  has  basis 

8^  *  »  1<1,3  i.p)  •  The  generalised  Bingham  matrix  distribution  (3.2) 

spans  Mj  •  **^,1,1  °*  <UjMnsion  v(p,p)  when  p^4  ,  and  «panc 

Kj  •  then  p  ■  3  .  a  basis  Is  providsd  by 

B2  ■  (x  ,  l<i<k<p,  l<.j<i,«p,  excluding  eases  corresponding  to  (3.3))  . 
when  ra  <  p,  corresponding  results  era  obtained  by  excluding  the  last  (p-r.) 
rows  of  X.  We  obtain  l*2  of  dimension  v(o,p)  with  basis 

(xljV,  lli£k£«,  113£*ip,  sxcluding  cssas 


(3.3)) 
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For  X-»hapee ,  and  the  distribution  (4.2)  us  procssd  similarly.  Me 
obtain  {x^x^.  l<_i£m',  1  <_  j£k  £p,  axcluding  casas  corresponding  to 
(4.2)),  of  dimansion  n(m',p).  For  hybrids,  and  tha  distribution  (5.2), 
a  basis  is  {x^,  a»1<i<m,  1<  j  <p}  •  t^j**** 

axcluding  casas  eorrasponding  to  (5.3)}  •  <  i<fc  <^m,  l<.3<i<.p» 

axcluding  casas  eorrasponding  to  (3.3)},  of  dimension  $ . 

Sines  Beran's  results  (1979)  apply  to  tha  canonical  exponential 
family  of  any  conpact  space,  his  formulae  (1.10)  for  estimators,  (5.5)  for 
approximate  tests,  and  their  extensions  to  Interval  estimation,  may  be  used 
in  large  samples  on  orientation  statistics  of  all  types,  provided  only 
that  suitable  robust  multivariate  density  estimates  are  available. 
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The  conventional  antlpodally  symmetric  Bingham  matrix  distribution  on  the 
Stlefel  manifold  Is  generalised.  Large  sample  maximum  likelihood  estima¬ 
tion  and  uniformity  tests  are  discussed,  and  a  parametric  model  for  axial 
orientations  (X*shapes)  Is  suggested.  A  generalisation  of  the  Khatrl-Nardla 
matrix  distribution  Is  developed  to  provide  a  model  suitable  for  hybrids 
(T-shapes).  Beran's  results  on  exponential  models  for  directional  data  are 
extended  to  orientation  statistics  to  provide  regression  estimators  and 
goodness-of-flt  tests  as  alternatives  to  maximum  likelihood  estimation  and 
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